塔里木河流域气候—水文过程的复杂性与非线性研究

【摘要】全球暖化现象正在世界各地以加速的节拍在发展,而伴随全球暖化的结果之一就是全球气候及水文循环正在急剧的改变,气候-水文过程因气温的上升及气候极端变化而有明显的调整。塔里木河流域作为一个干旱区典型,是气候水文变化的敏感区域,许多研究人员对其已经展开了大量的研究和讨论。对气候-水文过程的特征及其随时间演变的规律和内在动力学机制的研究,将有助于制定新的环保治理措施和发展更有效气候-水文预测预报方法,从而减缓气候-水文环境的恶化。气候-水文系统由多个子系统组成,并受到多圈层、多要素、多尺度的共同作用,子系统之间一种或多种方式发生非线性的内部或者外部的相互作用,导致它们不仅在时间上而且在空间上产生多种复杂形式的关联结构,形成了外有压迫、内有非线性耗散的复杂巨大系统,整体表现出复杂性。其影响因素有很多,例如,人类活动、地理位置、下垫面复杂的特征等等。但这些因素不是相互独立无关的,而是在各种时空尺度内产生复杂的非线性相互作用,从而使得塔里木河流域气候-水文过程的时间演化过程呈现出内在非线性外在复杂性的特征。针对上述问题,本文引入分形理论、混沌理论、复杂网络理论等并结合GIS技术对塔里木河流域的气候-水文过的复杂性与非线性特征进行了系统的定性与定量分析；深入研究了气候-水文因子某些复杂的时间演变过程和动力学机制；最后分析了气候-水文过程的临界自组织行为,对其存在的复杂性与非线性进行初步解释。本文主要的研究工作及其结论如下：1.首先应用重标极差法(R/S分析)、消除趋势波动分析法(DFA分析)以及谱分析法对日平均温度、日降水量、日平均相对湿度、日蒸发量和开都河日径流量序列进行分析,研究其长期持续性特征。结果表明,五个因子时间序列在较短的时间尺度上(约一年),具有明显的标度不变性,均表现出相类似的较高的长记忆性；而五个序列而在较长的时间尺度上(大于一年),不同因子时间尺度上的幂律特征存在着一些明显的差异。然后,借助消除趋势互相关性分析(DetrendedCross-CorrelationAnalysis,DCCA)揭示五个因子之间的长程互相性。结果发现,F2(S)与s的双对数曲线存在两个甚至三个无标度区间。最后在GIS空间分析技术支撑下,应用联合克利金插值法,分析研究了塔里木河流域气候-水文过程的长记忆强度的空间分布特征。2.根据混沌理论,对日平均温度、日降水量、日平均相对湿度、日蒸发量和开都河日径流量序列通过相空间重构,应用Rosenstein算法及G-P算法,分别计算了它们的序列的关联维数D、最大Lyapunov指数λ1以及Kolmogorov熵K等特征量。结果表明五个因子时间序列中存在明显的混沌特性,是非线性混沌动力系统演化的结果；塔里木河流域的气候-水文过程的变化同时具有确定性和随机性。最后在GIS空间分析技术支撑下,应用地统计理论,通过三个特征量分析了塔里木河流域的混沌特征和复杂度的分布。我们发现在气候变化复杂、具有径流、海拔较高的区域的值较大,混沌特性较强,复杂度较高,而受气候-水文过程影响较小的沙漠腹地区域的值较小,混沌特征较弱,复杂度较低,趋势明显。另外,通过关联维数D我们知道,D值越小(大),系统的层次越高(低),趋势越明显(不明显)。通过最大Lyapunov指数λ1,我们还知道五个因子的最大可预报尺度差异显著,低值地区最大可预报时间较长,高值区域较短；气候-水文因子系统的蝴蝶效应强弱程度表现为日平均相对湿度>日平均温度>日降水量≈日蒸发量,年径流量系统蝴蝶效应的平均强弱程度表现为和田河>叶尔羌河>阿克苏河≈开都河。3.应用两种多重分形方法(多重分形盒维数)对各气候-水文因子时间尺度和强度尺度上的多重分形结构进行分析。多重分形盒维数的分析表明,各因子在整个时间尺度上具有完全不同的多重分形特征,这些特征可以用三个具有明确物理意义的参数(B,Aa和△f)进行定量表征。进一步,本文将五个时间序列分为5个或者4个年代,分析了各年代因子的多重分形特征的变化,并结合塔里木河流域的气候水文变化进行了分析。最后在GIS空间分析技术支撑下,应用联合克利金插值法,分析研究了塔里木河流域气候-水文过程的多重分形强度的空间分布特征。4.根据粗粒化方法,将塔里木河流域的逐日平均温度、日降水量、日平均相对湿度、日蒸发量和开都河日径流量序列转化为由5个特征字符{R,r,e,d,D}构成的气候-水文因子符号序列。以符号序列中的125种3字串组成的温度、降水量、相对湿度、蒸发量、径流量的波动模态为网络的节点(即连续3d的因子波动组合),并按照时间顺序连边,构建有向加权的温度波动网络(TFN)、降水波动网络(PFN)、湿度波动网络(HFN)、蒸发波动网络(EFN)、径流波动网络(RFN)进而将五个因子的波动模态间的相互作用等综合信息蕴含于网络的拓扑结构之中。计算三种网络的度与度分布、聚群系数、最短平均路径长度等动力学统计量,从网络的角度对比研究五种序列内秉性质的差异,体现了五种气候-水文因子变化的复杂性。结果表明,五种因子波动网络均表现出很强的集聚性和较短的最小平均路径长度；TFN和RFN的度分布分别服从三段和双段幂律分布,是具有无标度特性的小世界网络；PFN、HFN和EFN兼具有无标度特性和小世界效应,既是无标度网络又是小世界网络。我们还发现在五种因子波动网络中,一些顶点的顶点度异常高,以这些顶点为代表的温度波动模态发生概率较大,例如,TFN中的节点RRR、dRR、ReR….。在TFN、HFN、EFN和RFN、PFN的重要节点中分别大都包含了R、r和r、e两种符号,这说明塔里木河流域的气候-水文过程的波动主要以上升为主,前三个因子比后两者上升的稍快。五个因子波动网络部分节点的介中向心性能力较强,4.5%(3.2%、3.2%、0.8%、3.2%)的节点承担了网络19.71%(19.71%、13.64%、3.4%、13.88%)的中介中心性功能,这些具有拓扑统计重要性的节点对于理解五个气候-水文因子波动的内在规律和波动信息的传递等有一定指导意义。根据5种顶点的平均度在塔里木河流域23个气象台站中的排行分布,将它们的主要作用区域依次定义为：一级度带、二级度带和三级度带。平均度、相似系数、最短平均路径长度和平均集群系数分布存在一定的区域特征,为通过复杂性来研究气候-水文变化的区域特征提供了理论基础。5.首先,五个气候-水文因子的强度-频度关系说明极端气候-水文事件与小因子值变化事件的过程均可归因于同一种动力学机制。然后,建立具有自身衰减因子的数值沙堆模型对气候-水文演变过程进行模拟。模型中仅仅通过改因子值的衰减系数,其结果就可以很好地解释五种气候-水文因子实际的强度-频度关系。模型研究的结果表明,因子自身的临界自组织行为是气候-水文过程时间演化的动力学根源,正是气候-水文因子自身的自组织临界性导致了其时间演变中存在幂律关系。
【作者】刘祖涵；
【导师】徐建华；
【作者基本信息】华东师范大学，地图学与地理信息系统，2014，博士
【关键词】复杂性；非线性；塔里木河流域；混沌；分形；复杂网络；自组织临界性；

【参考文献】
[1]刘利宝.全地面起重机传动方案的研究[D].太原科技大学,机械设计及理论,2014,硕士.
[2]林星锦.基于Android平台的无线订餐系统的设计与实现[D].北京交通大学,2013.
[3]田军军.催化金属表面结构和碳源对CVD法合成石墨烯的影响研究[D].重庆大学,化学工程与技术,2014,硕士.
[4]田立康.文本分类在电子取证领域的应用研究[D].北京化工大学,计算机应用技术,2013,硕士.
[5]欧景.双定子锥形永磁同步电机的弱磁研究[D].哈尔滨工业大学,电气工程,2013,硕士.
[6]仲伟来.农民家庭经营收入增长影响因素研究[D].扬州大学,产业经济学,2012,硕士.
[7]王新华.褐煤模型化合物与水相互作用及热解特性的理论研究[D].中国矿业大学,2014.
[8]李进.钢板剪力墙焊接残余应力及滞回性能研究[D].北京交通大学,2014.
[9]雪小峰.茂金属催化体系用于1-癸烯聚合制超高粘度PAO的研究[D].华东理工大学,油气储运工程,2014,硕士.
[10]崔文丽.SRS和RSS下加入辅助变量信息后总体均值的比率估计[D].华中师范大学,概率论与数理统计,2013,硕士.
[11]杨林.贲门周围血管离断术联合脾肾分流术治疗门静脉高压症疗效及Meta分析[D].第四军医大学,外科学（专业学位）,2014,硕士.
[12]尚德友,白玉山,张忠君,王起余.斜式旋耕的研究[J].农业机械学报.1992(03)
[13]王相林,王春芳.铁路客车检修信息系统的设计与应用[J].中国铁路.2002(07)
[14]董婵.种群动力学和信息安全中的数学模型研究[D].中北大学,应用数学,2014,硕士.
[15]金承焕,吴刚.实现业务创新的IT利器——TmaxSoftProductFactory[J].通信世界.2006(42)
[16]严海,王熙蕊,梁文博,孔令鹏.基于结构方程模型的通勤交通方式选择[J].北京工业大学学报,2015,04:590-596.
[17]李亚芳.基于技术维度的虚拟社区知识共享影响因素研究[D].华中师范大学,情报学,2013,硕士.
[18]李勋岳.建国前后职业报人的困境及其出路[D].华中科技大学,新闻学,2013,硕士.
[19]陈亚萍.灵武长枣果实糖代谢生理机制的研究[D].宁夏大学,植物学,2014,硕士.
[20]赵青玲.金属氧化物—石墨烯复合材料的合成及其电化学性能研究[D].湘潭大学,2014.
[21]彭翮翔.“国培计划”背景下中学英语教师培训中存在的问题研究[D].河北科技大学,英语语言文学,2012,硕士.
[22]葛祥旭.基于嵌入式系统的煤矿综采自动化通信控制系统设计与实现方法研究[D].山东大学,电子与通信工程（专业学位）,2013,硕士.
[23]张清林.汽车模型风洞关键影响因素研究及数据后处理软件开发[D].湖南大学,车辆工程,2012,硕士.
[24]肖军.悬索桥施工控制精细化分析及程序实现[D].重庆交通大学,桥梁与隧道工程,2012,硕士.
[25]安娜.中共八大至十二大全国党代会报告研究[D].中共中央党校,中共党史,2014,博士.
[26]曹方铖.江苏省高中化学优秀课堂教学案例的分析研究[D].扬州大学,学科教学,2011,硕士.
[27]周云.正畸牙齿移动过程中牙周膜成纤维细胞自噬作用的初步研究[D].第四军医大学,口腔正畸学（专业学位）,2014,硕士.
[28]何晴倩.美国选民“议题投票”行为的分析[D].华东师范大学,国际关系,2013,硕士.
[29]唐方萍.PAS-Na对锰致大鼠肝肾抗氧化酶和病理学改变的影响[D].广西医科大学,劳动卫生与环境卫生学,2013,硕士.
[30]张丽凤.虾味香精的开发及其在冷冻调理食品中的应用[D].广东海洋大学,水产品加工及贮藏工程,2013,硕士.
[31]李春宜.A宝钢冷轧板实现负公差轧制全自动化[J].冶金自动化,1996,02:10.
[32]樊平.跨文化视角下在华跨国乘用车研发中心的发展模式研究[D].上海外国语大学,工商管理,2012,硕士.
[33]蒋彦.以鲁棒性为目标的手势识别方法及其应用研究[D].济南大学,计算机软件与理论,2012,硕士.
[34]杨林华.农村扶贫模式创新研究[D].华中师范大学,公共管理,2014,硕士.
[35]孙彦宏.呼和浩特电信IP城域网的优化改造方案的设计[D].内蒙古大学,电子与通信工程,2014,硕士.
[36]滕志庆.本钢冷烧环冷机打滑检测保护装置的改造[J].冶金自动化,2013,01:73-75.
[37]王娟.四川工业遗产改造过程中工业元素景观再利用研究[D].西南交通大学,风景园林,2014,硕士.
[38]童思友,尚新民,石林光,王彦春.三维折射波静校正方法及应用[J].石油地球物理勘探,2002,S1:79-82+221.
[39]范天宝.复方黄芪对肉鸡免疫功能及抗氧化功能的影响[D].东北农业大学,临床兽医学,2013,硕士.
[40]孙光宇.大学生英语四级作文中模糊限制语失误的对比分析研究[D].吉林大学,外国语言学及应用语言学,2013,硕士.
[41]刘兴敏.西方绘画对中国现代工笔人物画的影响[D].东北师范大学,美术学,2012,硕士.
[42]贾慧丹.微纳流控芯片蛋白质富集分离新方法研究[D].华中科技大学,生物化学与分子生物学,2013,硕士.
[43]王启章.基于SaaS模式的基层医院合理用药监测与处方点评系统研究与应用[D].广东工业大学,工业工程,2013,硕士.
[44]顾芳芳.无限型Hopf π-代数的对偶[D].扬州大学,基础数学,2011,硕士.
[45]李艳.政府社会保障财政责任影响因素分析[D].南京财经大学,劳动经济学,2012,硕士.
[46]赵胜.基于B/S架构的Moodle网络教学平台的设计与实现[D].河北科技大学,计算机技术,2012,硕士.
[47]黄仁辉,张粒子,武亚光,庞博.中国电力金融市场的实现路径探讨[J].电力系统自动化,2010,11:54-60.
[48]郭翔.东江湖流域农田氮磷污染负荷特征研究[D].湖南科技大学,地图学与地理信息系统,2013,硕士.
[49]王秀丽.基于多层核心集的凝聚聚类算法[D].南京航空航天大学,应用数学,2012,硕士.
[50]张轶凡.货币政策对公司金融行为影响及其效果的研究[D].复旦大学,企业管理,2013,博士.